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一 种 基于 质量 评价 的 无 人 机 航拍 图 像 拼接 方法 
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摘 要 : 针对 大 视差 无 人 机 航拍 图 像 拼接 中 出 现 的 错位 、 伪 影 问题 ， 提 出 一 种 基于 质量 评价 寻找 最 佳 拼接 颖 的 方法 
QEB-U。 首先 通过 常规 拼接 缝 估计 得 到 初始 拼接 缝 ， 然 后 根据 无 人 机 航拍 图 像 的 特点 ， 综 合 考虑 结构 相似 性 、 色 差 、 
纹理 复杂 度 ， 提 出 一 种 新 的 质量 评价 函数 对 拼接 颖 上 的 每 个 像素 进行 评价 ， 进 而 根据 评价 结果 更 新 差异 代价 ， 之 后 
再 重新 估计 拼接 颖 。 估 计 和 评价 过 程 重复 执行 ,直到 拼接 颖 趋 于 稳定 则 停止 适 代 ,最 后 通过 梯度 融合 生成 最 终结 果 。 
实验 结果 表明 ， 所 提 方法 可 以 避免 大 视差 无 人 机 航拍 图 像 拼接 中 出 现 错位 、 伪 影 ， 优 于 目前 几 种 无 人 机 航拍 图 像 拼 
接 方法 ， 且 得 到 的 拼接 缝 优先 穿 过 道路 、 林 地 等 区 域 ， 更 加 符合 人 类 视觉 感知 ， 在 常见 的 图 像 清晰 度 评价 指标 上 均 
表现 良好 。 
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Quality evaluation-based image stitching for UAV images 


Wang Huan’, Zhou Yingyuea if, Wang Xueyuan’, Wang Xinyua 
(a. School of Information Engineering, b. Sichuan Provincial Key Laboratory of Robot Technology Used for Special 
Environment, Southwest University of Science & Technology, Mianyang Sichuan 621010, China) 


Abstract: Aiming at solving problems of dislocation and artifacts in the stitching of large parallax UAV images, this paper 
proposed a method to find the best stitching seam based on quality evaluation QEB-U. First, the initial seam was estimated 
by conventional seam estimation. Then, a new quality evaluation function was proposed for evaluating each pixel on the seam 
according to the characteristics of the UAV image like the structural similarity, color difference, and texture complexity. The 
difference cost was updated according to the evaluation result and the seam was acquired after that. The estimation and 
evaluation process was repeated until the seam becomes stable that the iteration was stopped. Finally, the stitching UAV image 
was generated by gradient fusion. The experimental results showed that the proposed method can avoid dislocation and 
artifacts in the stitching of large parallax UAV images, which was superior to several current UAV image stitching methods. 
Besides, the obtained stitching seam preferentially passed through roads, woodlands, etc. It was more in line with human 
visual perception and performed well in common image clarity evaluation indicators. 
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0 引言 作为 最 终结 果 。Lin 等 人 中 通过 自 适 应 特征 加 权 将 局 部 对 齐 计 
算 和 拼接 颖 估计 紧密 而 合 ， 并 通过 超 像素 特征 分 组 有 效 减 少 
随 着 无 人 机 遥感 技术 的 发 展 ， 其 研究 已 广泛 应 用 于 地 理 ”迭代 时 间 。Li 等 人 0 将 人 类 感知 的 非 线性 和 非 均匀 性 考虑 到 


岂 合 


言 恩 学 、 灾 害 监测 与 评估 、 国 防 等 领域 ， 尤 其 是 灾害 应 急 响 。 能 量 最 小 化 中 ， 用 显著 性 权重 来 模拟 人 了 眼 更 倾向 于 关注 的 对 
应 和 灾害 调查 中。 无 人 机 遥感 能 够 在 短 时 间 内 获得 目标 区 域 。 象 ， 使 得 拼接 结果 更 符合 人 类 视觉 感知 。Liao 等 人 0 指出 以 
最 直观 的 高 分 辨 率 图 像 ， 为 决策 提供 视觉 依据 。 然 而 受 无 人 ” 往 拼 接 颖 驱动 的 方法 都 是 评估 拼接 颖 的 平均 性 能 ， 无 法 找到 
机 飞行 高 度 和 相机 性 能 的 限制 ， 很 难 从 单 幅 无 人 机 航拍 图 像 ” 人 类 视觉 感知 上 最 佳 的 拼接 颖 ， 因 此 他 们 更 加 细致 的 评估 拼 
中 获取 整个 感 兴趣 区 域 的 景象 。 因 此 ， 有 必要 研究 将 具有 重 。 接 颖 上 的 每 一 个 像素 ， 并 引入 混合 成 本 来 评估 拼接 缝 沿线 的 
县 区 域 的 无 人 机 航拍 图 像 快速 拼接 成 具有 更 大 视野 的 高 分 辩 像素， 然后 根据 评估 结果 重新 估计 拼接 缝 ， 结 果 显 示 该 方法 
率 图 像 的 技术 外。 避免 了 伪 影 ， 并 且 总 体 上 效果 更 好 。Yuan 等 人 中 提出 一 种 基 
图 像 拼接 一 般 来 说 ， 有 两 种 方法 一 种 是 基于 单 应 变换 ”于 超 像素 的 无 人 机 航拍 图 像 拼接 方法 ， 拼 接 时 间 的 到 显著 提 
的 对 齐 技术 , 使 图 像 尽 可 能 精确 对 齐 E7; 另 一 种 就 是 基于 拼 ”，” 升 ， 并 在 能 量 函 数 中 融入 纹理 复杂 度 信息 ， 使 得 拼接 缝 穿 过 
接 颖 引导 的 切割 技术 ， 即 根据 预定 义 的 能 量 函 数 ， 寻 找 一 条 ”特定 的 低 纹 理 复杂 度 区 域 ,拼接 结果 更 加 符合 人 类 视觉 感知 。 
使 得 能 量 函 数 达 到 最 小 的 拼接 颖 Ps 。 基 于 单 应 变换 的 对 齐 Fang 等 人 (24 在 图 像 拼接 的 融合 阶段 进行 创新 , 提出 一 种 颜色 
技术 最 终 目 的 都 是 让 图 像 尽 可 能 对 齐 ， 但 是 当 遇 到 大 视差 的 混合 方法 来 消除 拼接 结果 中 的 颜色 不 一 致 问题 ， 并 利用 超 像 
环境 ， 该 类 方法 会 产生 局 部 伪 影 或 失真 的 情况 。 基 于 拼接 缝 ”。” 素 分 割 大 大 降低 计算 复杂 度 提高 拼接 效率 .Zhao 等 人 03 则 利 
引导 的 图 像 拼接 方式 最 早 是 由 Gao 等 人 提出 名 ， 首 先 通过 随 ”用 多 个 平面 代替 映射 点 进行 优化 ， 将 重 构 结 果 直 接 用 于 图 像 
机 抽样 一 致 性 算法 得 到 若干 单 应 变换 矩阵 ， 对 每 个 单 应 变换 融合， 实现 了 一 种 新 的 具有 较 高 精度 和 较 快 速度 的 航拍 图 像 
矩阵 计算 拼接 缝 ， 然 后 根据 质量 评价 结果 选择 最 好 的 拼接 缝 ”拼接 方法 。Pham 等 人 (首先 根据 航拍 图 像 的 地 理 位 置信 息 
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估计 出 重 辣 区域， 然后 将 特征 匹配 限制 在 该 区 域 ， 以 避免 不 
必要 的 计算 和 误差 累计 , 使 得 拼接 更 加 精确 快速 。 Ren 等 人 M1 
利用 航拍 图 像 的 位 置信 息 进行 粗 配 准 计算 出 全 局 变换 矩阵 ， 
并 通过 多 层 感知 机 将 全 局 变换 矩阵 进行 快速 校正 ， 因 此 该 方 
法 很 好 地 结合 了 基于 参数 和 基于 特征 的 方法 ， 实 现 了 拼接 
速度 和 精度 的 理想 折 中 。Guo 等 人 09 结 合 局 部 配 准 技术 提 
出 了 一 种 全 局 的 对 齐 策略 ， 使 得 拼接 结果 在 保持 形状 的 前 
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E(1)= >.D, 六 Spa ly, 
0W=2p()+ Bsrallohs) GO) 


其 中 N 代表 四 邻 域 ，4 是 2 的 相 邻 像素 点 。 常 规 拼接 颖 估计 
的 平滑 项 定义 为 

上 人 -ll (p)+L (gq)) (2) 

五 (EIAOmAO (3) 


提 下 获得 较 好 的 对 齐 精度 。Mo 等 人 07I 利 用 神经 网 络 进行 


中 差异 代价 代表 欧 几 里 德 色 差 d 中 。 数 据 项 D,(1,) 用 来 衡 


深度 特征 的 提取 与 匹配 ， 实 现 了 一 种 新 的 具有 和 鲁 棒 性 的 航 
拍 图 像 拼 接 方 法 ， 为 其 他 学 者 探索 结合 深度 学 习 的 图 像 拼 
帮助 。 
常 接 缝 驱动 的 方法 都 是 根据 预定 义 的 拼接 缝 质量 指标 
ee 但 是 这 些 方法 中 定义 的 质量 指标 都 
是 评估 拼接 颖 的 整体 质量 。 文献 [11] 已 证 明 , 在 许多 情况 下 , 整 
体质 量 最 好 的 拼接 颖 不 一 定 符合 人 类 视觉 感知 。 本 文 针对 大 视 
差 无 人 机 航拍 图 像 ， 提 出 一 种 新 的 拼接 颖 迭代 估计 方法 (quality 
evaluation-based image stitching for UAV ww QEB-U)， 其 
中 和 迭代 过 程 由 拼接 颖 上 的 每 个 像素 来 引导 。 首 先进 行 常规 # 
变 颖 估计 得 到 初始 拼接 缝 ， 然后 根据 无 人 机 航拍 图 像 的 特点 ， 
综合 考虑 结构 相似 性 、 色 差 、 纹 理 复 杂 度 的 特点 ， 提 出 一 种 
新 的 质量 评价 函数 对 拼接 颖 上 的 每 个 像素 进行 评价 ， 进 而 根 
评价 结果 更 新 差异 代价 并 重新 估计 拼接 颖 。 估 计 和 评价 过 
E 复 执行 ， 直 到 拼接 缝 趋 于 稳定 则 停止 迭代 ， 最 后 通过 梯 
kh 合 生成 最 终结 果 。 实 验 结果 表明 ， 所 提 方 法 对 大 视差 无 
L 航 拍 图 像 的 拼接 有 较 好 的 效果 ， 得 到 的 拼接 颖 优先 穿 过 
、 林 地 等 区 域 ， 更 加 符合 人 类 视觉 感知 ， 在 常见 的 图 像 
度 评 价 指标 上 均 表 现 良好 。 


所 提 方 法 


tH 


世 


融 几 


请 
E33 


为 了 找到 符合 人 类 视觉 感知 的 最 佳 拼 接 颖 ， 笔 者 更 加 细致 地 

去 评估 拼接 颖 上 的 每 个 像素 ， 提 出 一 种 新 的 基于 质量 评价 寻 

找 最 佳 拼 接 颖 的 方法 QEB-U, 在 解决 错位 伪 影 问题 的 同时 使 

ee eas 

1 .1 法 流程 
ee 


评估 
中 拼接 颖 上 的 上 
每 个 像素 


更 新 对 应 趋 于 稳定 
差异 代价 


最 佳 拼接 颖 


传统 
拼接 颖 估计 


图 1 QEB-U 方法 流程 图 
Fig.1 Flowchart of the QEB-U method 

拼接 颖 估计 外 得 到 初始 拼接 颖 ,之 后 对 
介 。 以 往 的 质量 评价 方法 都 是 基于 拼 
而 本 文 方法 则 是 对 拼接 颖 上 的 每 个 像 
素 进行 质 , 不 符合 期 待 的 像素 会 具有 较 大 成 本 , 符 
合 期 待 的 像素 则 县 有 较 小 成 本 。 成 本 较 大 的 像素 是 拼接 
颖 应 避免 穿 过 的 , 因此 后 续 拼 接 缝 重新 估计 中 , 在 能 量 函 


首先 采 
拼接 缝 进行 
接 颖 的 整 


J 波 
蔷 汪 
SS 


沁 型 


人 


-一 
法 


I 
Be 
下 


Se 
= 


~ 四 
坚 字 须 陨 天 


数 中 对 它们 进行 更 大 的 惩罚 。 当 拼接 缝 趋 于 稳定 则 停止 
和 迭代, 通过 梯度 融合 08 生 成 最 终结 果 。 下 面 将 关键 步骤 进 


行 详 细 说 明 。 
1.2 常规 拼接 缝 估 计 
设 两 幅 初 步 配 准 的 待 拼接 图 像 分 别 为 上 和 五 ，Zs(0D 是 


针对 大 视差 无 人 机 航拍 图 像 这 种 特殊 情况 下 的 图 像 拼接 ， 


量 


标签 分 配给 像素 2 的 代价 , 更 多 细节 请 参考 文献 [10]。 最 后 


通过 图 割 的 最 小 化 能 量 函 数 E(D) 来 估计 拼接 缝 。 


1.3 ”拼接 颖 迭代 估计 
1.3.1 质量 评价 


本 文 结合 无 人 机 航拍 图 像 特点 ， 综 合 考 虑 结构 相似 性 、 
色差 、 纹 理 复 杂 度 的 特点 ， 提 出 一 种 新 的 质量 评价 函数 。 首 


十 


讨 
先 重 登 区 域 未 对 齐 部 分 结构 上 通常 不 一 致 ， 因 此 在 质量 评价 
函数 中 引入 结构 相似 性 


指数 ( SSIM )P 来 评估 拼接 颖 。 其 次 引 


入 接 颖 两 侧 的 色差 来 更 加 细致 的 评估 拼接 颖 。 最 后 ， 对 于 无 


人 机 航拍 图 像 ， 发 现 符 合 人 类 视觉 感知 的 拼接 颖 更 倾向 穿 


域 纹 
氏 纹 
些 区 域 。 


自 


展示 了 
复杂 度 。 对 于 拼接 颖 上 的 像素 p; ， 混 合成 本 E(p;) 定义 式 为 


过 一 些 特定 区 域 ， 如 道路 、 林 地 、 湖 泊 等 。 根 据 这 些 区 
里 较 弱 的 事实 ， 在 质量 评价 中 引入 纹理 复杂 度 ， 使 得 
E 复 杂 区 域 的 像素 具有 更 小 的 成 本 ， 引 导 拼 接 颖 穿 过 
考虑 到 Gabor 滤波 器 5 可 以 很 好 地 表征 基础 纹理 
， 因 此 采用 Gabor 特征 描述 符 久 来 描述 纹理 复杂 度 。 图 
Gabor 特征 的 一 些 可 视 化 示例 ， 图 像 的 亮度 代表 纹理 


Ci 


kD JI 


E(pi)= 4 Boom (Pi): Epon (Pi) Evor (Pi) (4) 


中 4 是 一 个 常数 ， 保 持 ELP) 总 体 的 大 小 ，Bww Pi) 表示 以 


™ 


像素 pb; 为 中 心 的 局 部 


硬 片 在 了 4 与 1 之 间 的 SSIM 指数 ， 


Ewm(Pi) 表示 像素 bp 与 其 相 邻 像素 4; 的 色差 Eww (Pi) 表示 像 


素 p; 的 Gabor 特征 。 


SSIM 的 范围 为 [-1,1]， 当 两 个 局 部 区 域 越 相 似 ， SS 会 


趋向 于 1， 显 然 拼 接 颖 上 未 对 齐 的 像素 会 具有 较 大 的 值 。 
Eon (Pi) 的 表达 式 如 下 所 示 。 
Ena (p)= ee) (5) 
Ejom (Pi) We 
AA ET AA @ 


图 2 Gabor 特征 的 一 些 可 视 化 示例 


Some visualization examples of Gabor features 


Fig. 2 
1.3.2 更 新 差异 代价 


通常 ， 拼 接 颖 上 未 对 齐 的 像素 或 者 位 于 较 复杂 区 域 的 像 


一 个 标签 集 。 令 ? 代表 重 羞 区域 P 中 任意 一 个 像素 点 ， 拼 接 
颖 的 意义 是 将 一 个 标签 ,eIL 赋 给 每 个 像素 ? 。/,=0 代表? 的 
值 是 从 厂 复制 过 来 ， 而 =1 代 表 ? 的 值 是 从 了 复制 过 来 。 常 
规 拼接 颖 搜索 的 过 程 就 是 能 量 函 数 最 小 化 的 过 程 ， 能 量 函 数 


系 


有 较 大 的 混合 成 本 E(p;) 。 因 此 , 利用 像素 p; 的 混合 成 本 


E(p;) 对 差异 代价 14(p;) 进行 修正 ， 得 到 新 的 差异 代价 7(p): 


并 


PEeNV(sea) 


(7) 


otherwise 


天 、 [AGE La (pi) 
mo)-| Ti (pi) 


中 f(E(Pp)) 定 义 为 
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f(E(p))=e" Sm (8) 
其 中 N(seam) 是 包含 接 颖 在 内 的 带 状 区 域 ,通过 向 接 缝 左右 扩 
展 5 个 像素 生成 。o 是 一 个 控制 了 变化 幅度 的 常数 ，s 是 合 
理 成 本 的 闵 值 。 当 混合 成 本 E(p) 越 小 时 ，f(E(p)) 会 趋 近 于 
0， 新 的 差异 代价 元 (z) 也 会 越 小 ; 若 混 合成 本 E(p;) 越 大 时 ， 
JE(P)) 会 增 大 ， 新 的 差异 代价 元 (P) 也 会 越 大 。 然 后 利用 新 
的 差异 代价 重新 计算 能 量 函 数 , 并 通过 图 割 重新 估计 拼接 缝 。 
估计 和 评价 过 程 反复 进行 , 直到 拼接 缝 趋 于 稳定 ,这 里 的 “ 稳 
定 ” 意 味 着 当前 拼接 颖 可 以 完全 包含 在 先前 拼接 颖 的 带 状 
域 。 得 到 最 佳 拼接 缝 后 ， 通 过 梯度 融合 59 生成 最 终结 果 。 
1.4 算法 伪 代 码 
算法 1 QEB-U 

输入 : 两 幅 初 步 配 准 的 图 像 和 工 

输出 : 最终 的 拼接 颖 97 

1 初始 化 B= ; 

2 计算 差异 代价 矩阵 1 ,计算 平滑 项 5,4， 并 通 

数 E(1) 得 到 初始 拼接 颖 9 ; 

3 while Sig¢B 

4 计算 拼接 颖 % 上 每 个 像素 的 混合 成 本 E，; 

5 将 接 颖 8 扩展 为 带 状 区 域 N(5,)， 并 重新 计算 该 带 状 区 域 像素 的 

差异 代价 五 ; 

6 将 新 的 差异 代价 1 带 入 计算 平滑 项 sw ， 并 通过 图 割 最 小 化 能 量 

函数 E(1) 得 到 新 的 拼接 缝 9 ; 

7 令 B=BUN(S), $=5,; 

8 end 


人 


[| 


图 割 最 小 化 能 量 函 


说 


9 return S/;=5,; 
2 实验 


为 了 验证 本 文 所 提 方 法 的 有 效 性 ， 本 文选 用 公开 的 具有 
较 大 视差 的 无 人 机 航拍 图 像 来 进行 算法 性 能 测试 ， 其 中 数据 
来 源 于 文献 [2] 以 及 [24] 共 290 幅 图 像 。 实 验 中 ， 为 了 获得 最 
优 性 能 ， 将 相关 参数 设置 如 下 : 在 计算 局 部 片面 间 的 结构 相 
似 性 时 ， 

将 局 部 片面 大 小 设置 为 21x21， 式 (4) 中 的 4 设 为 10， 式 
(8) 中 的 c 和 = 设 为 5 和 0.12。 首 先 通过 文献 [9] 提 供 的 对 齐 方 
法 将 输入 图 像 初步 对 齐 ， 然 后 采用 所 提 方 法 进行 拼接 缝 欠 代 
估计 得 到 最 终 拼 接 终 ， 最 后 通过 梯度 融合 08 生 成 最 终结 果 。 
算法 实现 平台 为 Matlab R2020a， 主 机 为 2.9GHZ 八 核 CPU 
笔记 本 电脑 。 
2.1 与 其 他 拼接 方法 的 视觉 比较 

分 别 与 Zaragoza 等 提出 的 APAP(as-projective-as- 
Chang 等 提出 的 SPHP(shape-preserving half- 
projective)j51，Lin 等 提出 的 AANAP(Adaptive As-Natural-As- 
Possible)[d，Li 等 提出 的 ELA(parallax-tolerant image stitching 
based on robust elastic warping)25 以 及 Liao 等 提出 的 
SPW(single-perspective warps in natural image stitching)P9 方 法 
进行 比较 。 图 3(a) 是 一 组 具有 较 大 视差 的 机 场 无 人 机 航拍 图 
像 。 图 3(b) 是 APAP 算法 所 得 的 拼接 结果 , 图 3(c) 是 SPHP 算 
法 所 得 的 拼接 结果 ,图 3(d) 是 AANAP 算 法 所 得 的 拼接 结果 ， 
图 3(e) 是 ELA 算法 所 得 的 拼接 结果 ， 图 3(D 是 SPW 算法 所 
得 的 拼接 结果 ， 图 3(g) 是 所 提 算 法 得 到 的 拼接 结果 。 可 以 明 
显 看 出 图 3(b)-G9 的 机 标 和 机 尾 有 明显 的 伪 影 ( 红 框 部 分 )， 以 
及 明显 的 分 界 ( 绿 框 部 分 )， 而 本 文 所 提 方 法 获得 的 图 3(g) 所 
对 应 的 拼接 结果 清晰 无 伪 影 ， 且 更 加 自然 。 
图 4 是 一 组 具有 较 大 视差 的 房屋 无 人 机 航拍 图 像 及 几 种 
方法 的 拼接 结果 图 , 可 以 看 到 图 4(b)-( 了 的 房屋 出 现 了 伪 影 和 


possible)b] ， 


失真 模糊 ( 红 框 部 分 ) 图 4(b)(d)(e) 出 现 了 明显 分 界 ( 绿 框 部 分 )， 


而 图 4(g) 所 对 应 的 拼接 结果 清晰 自然 ， 无 伪 影 无 失真 。 
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像 


(c) SPHP (d)AANAP 


; 可 . Ce l ' \ L ' 
7 ~ FR ee eX 
EE : 
0 本 司 司 
(e)ELA (DSPW (g)QEB-U 
图 3 几 种 拼接 方法 对 于 有 较 大 视差 的 机 场 无 人 机 
航拍 图 像 的 拼接 结果 示意 图 


Fig.3 Schematic diagram of the stitching results several stitching 


methods for airport drone aerial images with large parallax 


Ge (OSPHP ANA 


Ko FoR eed ey 
(e)ELA (DSPW (g)QEB-U 
图 4 几 种 拼接 方法 对 于 有 较 大 视差 的 房屋 无 人 机 

航拍 图 像 的 拼接 结果 示意 图 


Fig.4 Schematic diagram of the stitching results Several stitching 


methods for aerial images of house with large parallax 
2.2 与 其 他 拼接 缝 搜索 算法 的 视觉 比较 
为 了 验证 本 文 所 提 拼 接 颖 搜索 策略 的 有 效 性 ， 本 文 与 最 
新 的 拼接 颖 搜索 算法 进行 了 对 比 ， 主 要 包括 Perception- 
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based(PB)010 ，Qnuality evaluation-based(QEB)00 和 Fast and ”都 在 一 分 钟 以 上 。 本 文 所 提 方 法 在 拼接 效率 上 与 PB 和 QEB 

Tobust(FARSE)P7 方 法 。 相似 ， 但 不 及 FARSE， 因 此 所 提 方 法 的 拼接 效率 还 有 待 提 高 ， 后 
如 图 5 所 示 ， 红 线 标记 了 各 个 方法 检测 到 的 最 佳 拼 接 缝 续 考虑 使 用 超 像 素 分 割 | 来 减少 迭代 计算 量 ， 从 而 提高 拼接 效率 。 

位 置 。 可 以 看 到 ， 由 PB, QEB, FARSE 方法 寻找 到 的 最 佳 拼 ”2.4 客观 指标 评价 


接 颖 会 穿 过 突出 物体 , 例如 房屋 等 建筑 物 ( 绿 框 部 分 ), 而 所 提 为 了 更 加 精准 客观 地 比较 不 同方 法 拼接 后 的 图 像 质量 优 
方法 寻找 到 的 拼接 缝 可 以 避 开 房屋 等 突出 物体 ， 优 先 穿 过 道 劣 , 本 文 使 用 一 些 常见 的 图 像 清 晰 度 评价 指标 , 如 Brenner 梯 
路 、 林 地 等 ， 使 得 拼接 效果 更 加 符合 人 类 视觉 感知 。 度 函 数 、Tenengrad 梯度 函数 、Laplacian 梯度 函数 、 灰 度 方差 
2.3 与 其 他 拼接 缝 搜索 算法 的 效率 比较 函数 以 及 能 量 梯度 函数 。 上 述 几 种 函数 计算 的 值 越 大 ， 表 明 

同时 ， 本 文 与 最 新 的 拼接 颖 搜索 算法 PB、QEB、FARSE ”拼接 图 像 的 清晰 度 越 高 ， 拼 接 效 果 越 好 。 下 面 对 这 几 个 函数 


在 拼接 时 间 上 进行 了 对 比 ， 本 文选 用 [24] 提 供 的 三 组 分 辨 率 ”的 计算 做 简要 描述 。 


为 3680*2456 的 高 清 航 拍 图 像 进行 对 比 实验 分 析 。 表 1 是 三 Brenner 梯度 函数 Ds(f) 定义 为 
组 实验 的 具体 数据 ， 可 以 看 出 在 对 和 高 分 辩 率 航拍 图 像 进 Ds =D, Dt -fF (9) 


行 拼接 时 , PB、QEB、FARSE 以 及 本 文 所 提 方 法 的 拼接 用 时 其 中 f (wy) 表示 图 像 1 对 应 像素 点 (%y) 的 灰 度 值 。 


(b)QEB (CJFARSE (d)QEB-U 
图 5 不 同 拼接 颖 搜索 算法 的 最 佳 拼 接 颖 位 置 比较 


Fig.5 Comparison of optimal seam positions for different seam search algorithms 


表 1 不 同 拼接 颖 搜索 算法 的 拼接 效率 对 比 Tenengrad 梯度 函数 D:(f) 定义 如 下 : 
Tab. 1 Comparison of stitching efficiency of DP=DD EY Ty)>T) (10) 
different stitching seam search algorithms T(x,y) 的 定义 如 下 : 
分 组 不 同 拼接 算法 ”拼接 用 时 (s) TC) = NO D+)) (11) 
ge PB 256 其 中 了 是 给 定 的 边缘 检测 阔 值 ， G6 和 9， 分 别 是 像素 点 (x,y) 
Ws QEB 162 处 Sobel 水 平和 垂直 方向 边缘 检测 算 了 的 郑 职 。 
8 FARSE 65 Laplacian 梯度 函数 D.(f) 定义 如 下 : 
Oe 680 2945p QEB-U 158 DOD= 了 ,Lo (Ly)>7) (12) 
晤 前 PB 146 其 中 L(x,y) 是 像素 点 (%y) 处 Laplacian 算 子 的 卷 积 。 
图 像 大 小 ，472MB 与 4.62MB QEB 300 灰 度 方差 函数 D;(f) 定义 如 下 : 
Hy FARSE 76 |f G3) -fy -DI+ 
分 辨 率 :3680*2456 oe 5 wi ] (13) 
第 三 组 同 148 能 量 梯度 函数 D:(7) 定义 如 下 : 
图 像 大 小 : 3.75MB 与 3.73MB SE Ee [f(x+l,y)-f (x,y) + 
分 辨 率 ，3680*2456 Ee p05 (wy+D)-f (wy) | "9 
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对 图 3 和 4 的 拼接 结果 分 别 进 行 清晰 度 指 标 计算 并 进行 
线性 归 一 化 处 理 ， 实 验 结果 如 表 2、3 所 示 。 可 以 看 出 : 本 文 
法 几乎 在 所 有 评价 指标 上 都 比 其 他 方法 具有 更 高 的 数值 ， 

这 意味 着 所 提出 的 方法 可 以 获得 更 加 清晰 的 拼接 图 像 。 


表 2 图 3 的 清晰 度 测试 结果 比较 
Tab. 2 Comparison of clarity test results in figure 3 


不 同 拼接 “Brenner Tenengrad Laplacian 灰 度 能 量 
算法 梯度 函数 ”梯度 函数 ”梯度 函数 ”方差 ”梯度 
APAP 0.7660 0.5467 0.9291 0.8143 0 
SPHP 0.4057 1 0 0 1 
AANAP 0.7912 0.7606 1 1 0.0762 
ELA 0 0.8999 0.4348 0.4596 0.7053 
SPW 0.0785 0 0.4446 0.1534 0.1071 
QEB-U 1 0.9536 0.8148 0.8519 0.2018 


表 3 图 4 的 清晰 度 测试 结果 比较 
Tab.3 Comparison of clarity test results in figure 4 


不 同 拼接 ”Brenner ”Tenengrad Laplacian 灰 度 能 量 
算法 梯度 函数 ”梯度 函数 ”梯度 函数 ”方差 梯度 
APAP 0.7148 0.6929 0.7668 0.7299 0.7025 
SPHP 0 0 0 0 0 
AANAP 0.0490 0.0530 0.0256 0.0396 0.0334 
ELA 0.5079 0.5936 0.2595 0.4095 1 
SPW 0.7052 0.7285 0.6466 0.6828 0.9182 
QEB-U 1 1 1 1 0.9415 
3 ”结束 语 
本 文 针 对 无 人 机 航拍 图 像 拼 接 提出 一 种 基于 质量 评价 寻 
找 图 像 


最 佳 拼接 颖 的 策略 ， 不 仅 可 以 解决 大 视差 无 人 机 航拍 医 


人 
提 


觉 


接 中 错位 伪 影 的 问题 ， 而 且 得 到 的 拼接 缝 在 内 容 复 杂 自 
机 航拍 图 像 中 优先 穿 过 道路 、 林 地 等 区 域 。 实 验 表 明 ， 所 
方法 适合 于 无 人 机 航拍 图 像 拼接 ， 拼 接 结果 更 符合 人 类 视 
感知 。 
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